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Abstract 

During  the  course  of  this  grant  improved  algorithms  for  simulating  photoelectron  spectra 
were  developed  and  used  to  simulate  the  photoelectron  spectra  of  the  cyclopentadienyl  (CH)5 
substitutional  isomers  triazolyl  (CH)2N3  and  pyrrolyl  (CH)4N.  The  programmatic  significance  of 
this  task  lies  in  the  importance  of  the  parent  azoles,  (CH)„N5nH,  and  their  deprotonated 
counterparts,  the  azolides,  as  fundamental  units  in  ionic  liquids,  explosives  and  fuels.  Significant 
progress  was  also  made  in  a  study  of  the  nonadiabatic  photodissociation  of  furazan,  (CH)2N20 
another  five  member  ring  species.  Furazan-based,  high  nitrogen  content  energetic  materials 
have  recently  received  significant  attention  for  applications  as  explosives,  fuels,  and  propellants 
due  to  their  good  heat  resistivity,  high  heat  of  formation,  low  sensitivity,  and  good  detonation 
performance.  For  these  furazan  based  materials  photo-induced  dissociation  of  the  electronically 
excited  species  is  relevant  to  their  denotation. 
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Narrative 


I.  Overview 

The  principal  task  we  have  addressed  in  the  course  of  this  grant  is  the  determination  of 
the  negative  ion  photoelectron  spectra  of  azolides,  precursor  anions  to  the  azolyls,  nitrogen 
heterocycles  of  the  form  (CH)5_„Nn,  1  <  n  <  5,  whose  vibronic  spectra  are  observed  in  the 
experiment,  azolide  +  hv  — >  azolyl  +  e  .  The  programmatic  significance  of  this  task  lies  in  the 
importance  of  the  parent  azoles,  (CH)„N5_nH,  and  their  deprotonated  counterparts,  the  azolides,  as 
fundamental  units  in  ionic  liquids,1,2  explosives3  and  fuels.4  The  thermochemistry  of  these  azoles 
and  negative  ion  photoelectron  spectra  of  the  azolides,  have  been  reported  by  Lineberger's  group 
at  the  University  of  Colorado.5'8  The  photoelectron  spectra  are  quite  complicated  and  most  are 
poorly  described  by  adiabatic  techniques.6,7'9  In  order  to  carry  out  the  simulations  it  was 
necessary  to  make  significant  upgrades10,11  to  the  algorithms  we  had  only  recently  developed12'14 
to  describe  nonadiabatic  photoelectron  spectra.  The  complexity  is  attributable  to  the  existence  of 
low-lying  conical  intersections  of  two  and  three  electronic  states  which  are  not  required  by 
symmetry.7,15,16  The  existence  of  low-lying  accidental  three  state  conical  intersections  in  these 
molecules,  which  was  completely  unprecedented,  was  established  in  the  course  of  earlier 
AFOSR  funded  work.15  The  computational  advances,  and  the  computed  spectra  and  their 
significance  are  described  below. 

Nonadiabatic  photodissociation  can  play  a  role  in  detonation  of  energetic  materials  when 
the  ignition  processes  involves  sparks,  shocks  or  laser  heating.17  We  considered,  in  a  project 
which  is  ongoing,  the  nonadiabatic  photodissociation  of  furazan,  C2H2N20*  -»  NO  +  C2H2N. 
Our  work  brings  new  insights  into  a  problem  carefully  studied  experimentally  and 
computationally  in  the  laboratory  of  Bernstein18  (Colorado  State). 

II.  Improved  Algorithms  for  KDC  Theory 

A.  An  Enhanced  Quadratic  Vibronic  Coupling  Model-Existing  Tools 

The  spectra  computed  during  the  course  of  this  grant,  those  of  pyrrolide  and  1,2,3- 
triazolide,  were  to  be  computed  using  our  second  order  implementation  of  the  multimode 
vibronic  model  of  Koppel,  Domcke  and  Cederbaum(KDC  Theory), 19,20  introduced  30  years  ago 
and  extended  and  improved  since  then  by  both  those  workers21'24  and  others.7,12"14,25  In  KDC 
theory,  the  computation  of  the  photoelectron  spectrum  proceeds  in  two  steps.  In  this  presentation 


we  use  the  language  of  negative  ion  photodetachment,  although  the  approach  is  also  suitable  for 
photoionization.  In  the  first  step  the  Nsuae  coupled  adiabatic  electronic  states  of  the  neutral  are 
described  using  a  coupled  quasi-diabatic  state  Hamiltonian,  Hd,  constructed  in  N"n  nonredundant 
internal  coordinates,  Q 

Hfj( Q)  =<  ,Q)  I  I  WrUr""'  ,Q)  >=  £,*( Q°)<5„  +  '%V™JQk  +1/2 '2,V™JQkQm  (1) 

k  k,m 

where  1  <  i  <  j  <  Wtate,  W.  ’d  is  a  quasi-diabatic  electronic  state,  the  coordinates  of  the  If1 
electrons  are  denoted  rN  =  (iq,  r2,  ...  .  r  e()  and  E]  is  the  adiabatic  energy  of  the  fh  state.  Note 

that  Hd  contains  terms  through  second  order  in  Q,  and  the  resulting  KDC  method  is  denoted  the 
quadratic  vibronic  coupling  (QVC)  model.  In  our  extension  of  the  KDC-QVC  method,  m.4  is 


not  explicitly  determined  but  rather  the  V(l)  are  determined  by  requiring  the  adiabatic  energies, 
energy  gradients  and  derivative  couplings  obtained  from  Hd  to  reproduce,  as  well  as  possible,  the 


1 


corresponding  ab  initio  data.  The  unique  aspect  of  our  approach  is  that  derivative  coupling 

information  is  used  in  the  construction  of  Hd.  This  results  in  a  representation  which  is 

1  ^ 

quantifiably  quasi-diabatic  in  a  least  squares  sense. 

In  the  second  step,  the  total  vibronic  Schrodinger  equation 

(Ht-E])WTj  =0  (2a) 

where  HT=  Twc+  HeN  ,  and  Tnuc  is  the  nuclear  kinetic  energy  operator,  is  solved.  The  eigenstates 
are  expanded  in  a  vibronic  basis 

N  state 

^(rw-,Q)=  2^(r^';Q)2^  ln>  (2b) 

(=1  n 

where  the  vibrational  basis,  |n>,  is  given  in  the  multimode  representation  by24 

Nmt 

|n(Q))-n*»,Ga)  for  0  <  /?,  <  Mi  (2c) 

i=  1 

and  x'n  is  the  nth  harmonic  oscillator  function  for  the  ith  mode  of  the  neutral.  In  this  basis  eq.  (2a) 
becomes: 

[Hr  -  £^l]c*  =  0  (3a) 

where  «,'«>-(»  I  T""StJ  +  I  m)Q  .  (3b) 

The  multimode  vibronic  basis  has  dimension 

Nmt 

Nmm=  NsmteNvib  where  Nvib  =  \\Mt  (4) 

i=  1 


which  can  easily  exceed  a  billion  terms,  depending  on  the  choice  of  Mt.  However  given  the 
polynomial  form  of  Hd  and  harmonic  oscillator  products  in  n(Q),  HT  has  structured  sparsity.24 
The  number  of  nonzero  elements  grows  only  linearly14  with  NMM  and  the  location  of  the  zeros  in 
HT  is  readily  identified.  cK  is  determined  from  eq.  (3a)  using  a  Lanczos  based 
“diagonalization.”26  As  part  of  an  earlier  AFOSR  grant,  we  developed14  an  open-ended  fine¬ 
grained  parallel  Lanczos  algorithm,  which  partitions  the  NMM  dimensional  vector,  cK  among  IVproc 
processors.  This  both  reduces  the  time  to  solution  of  eq.  (3a)  and  increases  the  maximum  NMM 
which  is  practical.  As  currently  implemented  we  can  handle  NMM~  1010,  which  is  one  to  two 
orders  of  magnitude  larger  than  that  possible  with  single  processor  codes. 

The  simulated  photoelectron  spectrum  is  given  by  the  spectral  intensity  distribution 
function,  P(E)  where 

P(E)  =2ji2\Bk\2S(E-(Etk-E0))  where  BK=(fi^\wTK)  .  (5) 

K 


Here,  /t  is  a  many  electron  dipole  moment  operator,^  is  the  vibronic  wave  function  of  the 
negative  ion  which  is  assumed  to  be  in  its  ground  vibrational  level  and  well  described  in  the 
adiabatic  state-harmonic  oscillator  (the  quadratic,  adiabatic)  approximation,  that  is 


Ar+i 


,Q)  =  ip-'"(rw<',+l;Q) 1 0  >=  W-«(r"''+l;Q)24n  I  n  > 


(6a) 


where  ^"(r^Q)  is  the  ground  adiabatic  electronic  state  of  the  anion,  |o)  denotes  the 
ground  vibrational  level,  the  overbar  indicates  that  the  vibrational  coordinates  and  harmonic 
oscillator  functions,  are  appropriate  for  the  anion  or  precursor  molecule  and  may  differ  from 
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those  of  the  neutral.  The  are  multidimensional  Franck-Condon  factors 

In""  n""  \ 

C  =<  0  I  m  >  where  (n  |  m)  =  ( II E,.  (ft )  II  xL,  0 G,  )j  •  (6b) 

The  spectral  intensity,  \BK\2  for  the  transition  from  the  ground  level  of  the  negative  ion,  given  by 
eq,  (6a),  to  the  K1"  state  of  the  neutral,  given  by  eqs.  (2b,  2c)  is  readily  shown  to  be:12 


bk  =  ~  Ho 


do/'oj,  )c/n 


where 


=(md0..;n)Q 


do,v(Q)  =  (wr"(rJve'+1;Q)|/)|  2ip/’"(r^  ;Q)>  •  (8b) 

/  r 

The  c*,  are  obtained  from  eq.  (3a)  using  the  open-ended  Lanczos  algorithm.  The  second 
equality  in  eq.  (7)  or  (8a)  is  valid  when  the  Q  dependence  of  u0  ;(Q)  can  be  neglected.  BK  can  be 
determined  without  explicit  determination  of  cK]a  if  the  seed  vector  in  the  Lanczos  procedure  is 
properly  chosen.712,24 

We  also  showed12  that  with  a  judicious  choice  of  the  x'm  (ecl-  2c),  NMM  can  be  reduced 
considerably,  by  factors  of  30  to  300  in  our  test  case,12  further  extending  the  range  of  problems 
that  can  be  treated.  In  order  to  exploit  this  basis  set  flexibility  the  Franck-Condon  factors  in  eq. 
(6b)  are  required.  We  developed12  an  algorithm  which  uses  fine  grained  parallelism  to  efficiently 
and  exactly  determine  the  millions  to  billions  of  multidimensional  Franck-Condon  overlap 
integrals.19,27  31 . 

The  electronic  transition  dipole  moment  in  eq.  (8b)  deserves  further  comment.  Note  that, 
eq.  (8b)  is  significantly  oversimplified,  since  the  anion  wave  function  contains  Nel+ 1  electrons 
while  that  of  the  neutral  contains  only  Nel electrons.  The  discrepancy  is  attributable  to  the  failure 
to  incorporate  a  description  the  detached  electron.  This  issue  is  addressed  in  the  research 
proposal  submitted  as  the  renewal  of  this  grant. 

B.  Higher  order  contributions  to  Hd 

As  discussed  below,  the  QVC  model  described  above,  currently  the  workhorse  of  the 
field,  worked  quite  well  for  the  photoelectron  spectrum  of  pyrrolide  [  (CH)4N  ]  but  proved 
inadequate  for  the  description  of  the  spectrum  of  triazolide  [  (CH)2N3  ].  Consequently  during 
the  current  performance  period,  the  domain  of  utility  and  the  accuracy  of  were  extended  by 
including  higher  order  terms.11  Our  Hd,  through  fourth  order,  is  given  by 

+ 2'/.<wa  -  «T’<Q)  <9a> 


«5“’<Q)  =  4  2v‘^Q‘S'Q-+7i  aae.  w 

*  k,l,m= 1  *  j,k,l,m= 1 

The  algorithm,  to  determine  this  Hd  has  been  discussed  in  detail  in  Ref.  11 .  Unlike  other 
implementations  of  higher  order  terms,32’33  our  implementation,  does  not  require  the 
simplifications  due  to  high  symmetry.  As  with  the  Hd  in  eq.  (1)  this  representation  is 
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quantifiably  quasi-diabatic  in  a  least  squares  sense.11 

Our  simulation  of  the  triazolide  spectrum  described  below10  which  required  V  through  4th 
order,  provided  insights  into  the  need  for  higher  order  terms.  The  fourth  order  terms  were 
required,  not  because  Ns’“,e  =  4  was  needed,  for  the  other  azolyls  treated  Ns,a,e  =  2  (Refs. 16,34)  or  3 
(Refs.  7,3S)  had  been  adequate  but  rather  because  of  the  comparatively  large  geometrical 
separation  of  key  features  of  the  adiabatic  potential  energy  surfaces,  minima,  transition  states  and 
conical  intersection  seams.11 

C.  A  general  open-ended  Lanczos  procedure  for  Hd  with  Higher  Order  Terms. 

The  polynomial  form  in  eqs.  (9a, 9b)  preserves  the  structured  sparsity  of  HT  in  eq.  (3a). 
As  part  of  our  AFOSR  funded  research,10  we  used  this  property  to  develop  a  fine-grained  parallel 
version  of  the  Lanczos  "diagonalization"  algorithm  that  can  handle  arbitrary  V(n)  contributions. 
This  algorithm,  which  is  now  in  routine  use,  allows  us  to  exploit  the  enhanced  capability  of  the 
higher  order  Hd  to  achieve  accurate  simulations  in  computationally  challenging  systems. 

III.  Azolide  Photoelectron  Spectra. 

We  now  turn  to  the  results  of  the  simulations  of  the  azolide  photoelectron  spectra  which 
reveal  the  spectra  of  the  azolyls.  The  (CH)mN5.m“  heterocycle  anions  are  of  programmatic 
interest  as  potential  anionic  components  of  ionic  liquids  and  as  tools  to  infer  properties  of 
pentaazolide,  N<f.9,36  The  thermodynamics  of  these  anions  was  studied  both  experimentally  and 
computationally  by  Lineberger's  group  at  the  University  of  Colorado.5-8,37,38 

The  azolyls  are  substitutional  isomers  of  the  D5h  cyclodentadienyl  radical,  (CH)5  whose 
2Ej"  ground  state  is  Jahn-Teller  distorted.39-41  Replacing  CH  groups  with  N  atoms  both  reduces 
the  point  group  to  C2v  and  introduces  in-plane  lone  pair  orbitals,  excitations  from  which  lead  to 
low-lying  excited  states.  The  2E,"  symmetry-required  conical  intersection  seam  in 
cyclopentadienyl  becomes  a  symmetry-allowed  seam  in  an  azolyl  and  the  additional  low-lying 
excited  states  can  lead  to  additional  conical  intersection  seams.  Since  our  enhanced  QVC 
algorithms,  discussed  in  Section  II A,  had  previously  enabled  us  to  determine  the  photoelectron 
spectrum  of  the  pyrazolide  -  pyrazolyl  system  [(CH),N2J  (Ref. 35)  which  had,  at  that  time,  been 
declared  'noncomputable'  using,  then  state  of  the  art  time-independent  QVC  techniques,7  we 
expected  the  nonadiabatic  pyrrolide  and  triazolide  simulations  to  be  readily  amenable  to  our 
QVC  treatment.  This  did  not  turn  out  to  be  the  case. 


A.  Photoelectron  spectrum  of  the  pyrrolide-pyrrolyl  system: 


The  most  acidic  site  on  pyrrole,  C4NH5  is  the  nitrogen.9  Deprotonation  at  this  site 
produces  the  anion  pyrrolide,  (CH4)N“,  whose  photoelectron  spectrum  reveals  the  vibronic 
spectrum  of  pyrrolyl.  In  pyrrolyl,  pictured  above,  the  2E1"  ground  state  of  cyclopentadienyl  is 
split  into  a  A2  ground  state  and  a  Bi  excited  state.  The  symmetry-required  Ei"  conical 
intersection  seam  in  cyclopentadienyl  becomes  a  symmetry-allowed  2A2-  2Bi  accidental  seam  of 
conical  intersection,  whose  minimum  energy  crossing  point  is  4436  cm1  above  the  2A2  ground 
state  minimum.34  The  minimum  energy  C2v  nuclear  configuration  where  the  wave  function 
carries  a  2Bi  irreducible  representation  is  a  saddle  point.  Fig.l  presents  the  pyrrolide 
photoelectron  spectra  determined  experimentally9  and  simulated34  using  our  nonadiabatic  two 
state  QVC  method.  Domcke  and  coworkers16  have  also  reported  a  nonadiabatic  simulation. 
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Figure  1:  Pyrrolide  photoelectron  spectrum.  (a)Experimentally  measured9  (b)Computed  nonadiabatic  simulation34 
From  this  figure  is  seen  that  when  the  electron  binding  energy  (eBE)  is  less  than  -  2.4  eV  the 
agreement  between  theory  and  experiment  is  quite  good.  Similar  conclusions  hold  for  a  purely 
adiabatic  determination  of  the  photoelectron  spectrum  based  on  the  A2  state.  However,  as 
originally  noted  in  Ref.  9,  attempts  to  include  the  effects  of  the  2Bi  state  using  an  adiabatic 
approach  produce  "a  prominent,  well-resolved  photoelectron  spectrum  beginning  about  0.5  eV 
above  the  ground-state  origin."  But  "no  such  spectrum  is  observed."9  As  discussed  in  Refs.  9 
and  16  and  confirmed  by  our  calculations,34  this  a  consequence  of  the  nonadiabatic  interactions 
induced  by  the  2A2  -  2Bi  seam  of  conical  intersection.  This  observation  is  not  entirely  surprising 
since,  as  noted  above,  the  minimum  energy  conical  intersection  is  -0.55  eV  (4436  cm'1)  above 
the  ground  state  origin. 

The  good  agreement  of  the  intensities,  peak  heights,  between  the  measured  and  our 
simulated  spectrum  for  electron  binding  energy  (eBE)  >  ~  2.6  eV  is  a  consequence  of  two 
factors:  the  inclusion  of  nonadiabatic  interactions  through  Hd  and;  adjusting  r  =  u  2 ,  /  u  2 u 

(J ,  ^2  U ,  lj  y 

24 

from  the  standard  choice  r  =  1,  to  r  =  2.  It  was  need  for  this  empirical  correction  to  the 
transition  dipole  moments  that  served  to  motivate  the  approximate  first  principles  determination 
of  , u0  j  discussed  in  our  AFOSR  renewal  proposal. 

A  question  of  fundamental  importance,  which  has  received  only  limited  attention,42’43  is 
how  far  below  the  minimum  energy  crossing,  are  the  effects  of  the  conical  intersection  seam  felt. 
Analysis  of  the  vicinity  of  the  peak  labeled  B  in  Fig.  1  at  eBE  2.257  eV  (905  cm"1)  evinces  a 
shoulder  at  2.229  eV,  only  675  cm'1,  0.084  eY,  above  threshold  and  -3760  cm"1  below  the 
minimum  energy  crossing.  A  symmetry  analysis  of  our  wave  functions  shows  that  this  shoulder, 
which  is  absent  in  the  adiabatic  simulation,  is  due  to  nonadiabatic  coupling  to  the  Bi  diabat. 
With  the  advent  of  high  resolution  SEVI  and  PFI-PE  spectroscopy  such  low  intensity  peaks 
should  be  readily  discerned.  Given  that  the  feature  in  question  is  -3760  cm'1  below  the 
minimum  energy  crossing  its  nonadiabatic  character  is  quite  surprising.  The  extent  of 
nonadiabatic  intrusion  into  ostensibly  adiabatic  regions  of  photoelectron  spectra  is  an  important 
fundamental  issue  of  practical  import.  This  question  will  be  addressed  in  our  AFOSR  renewal 
proposal. 

B.  Photoelectron  spectrum  of  the  triazolide-triazolyl  system: 


The  1,2,3-triazolide  anion,  (CH)2N3  ,  is  obtained  by  deprotonaton  of  1H- 1,2 ,3 -triazole  at 
the  most  acidic  site,  the  nitrogen  atom.  This  system  has  been  studied  experimentally  and 
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computationally  in  Lineberger's  laboratory.38  Lineberger  observed  that  the  experimental 
spectrum  could  not  be  simulated  using  a  quadratic,  adiabatic  description  of  the  neutral.38  The 
failure  of  the  quadratic,  adiabatic  description  of  1,2,3-triazolyl  was  surprising  in  view  of  the  fact 
that  the  minimum  energy  point  on  the  conical  intersection  seam  is  3737  cm'1  (0.463  eV)  above 
the  ground  state  minimum.10  In  the  case  of  pyrrolide  we  were  able  to  identify  small  nonadiabatic 
effects  in  the  experimental  spectrum  at  energies  that  extend  well  below  the  minimum  energy 
point  on  the  conical  intersection  seam.  Attempts  to  perform  a  similar  analysis  for  1,2,3- 
triazolide  were  complicated  by  the  fact  that  the  nonadiabatic  effects  lie  in  a  spectral  region  where 
hot  band  sequences  exist.  Thus  the  simulation  of  photoelectron  spectrum  of  1 ,2,3-  triazolide  was 
particularly  challenging  since  it  required  us  to  address  a  combination  of  issues. 

Fig.  2  compares  our  nonadiabatic  simulation,  with  the  experimentally  measured 
spectrum.  The  agreement  is,  as  discussed  in  Ref. 10,  quite  reasonable.  Below  we  discuss  how,  in 
general,  this  agreement  was  achieved,  and  the  spectrum  near  peak  A,  the  origin  peak  identified 
by  Lineberger. 


eBE(eV) 


(a)  (b) 

Figure  2:  Triazolide  photoelectron  spectrum,  (a)  Experimental  spectrum.38  (b).  Simulation  of  based  on  Hd 
including  Hd  <3,4>.  Contributions  originating  from  five  vibrational  levels  of  the  anion  are  included:  ground  state(red)  , 
remainder  have  one  phonon  in  0>7  (purple),  CP8(green),  fO ,/blue),  and  CO  ,,/yellow),  respectively. 


The  simulations  of  the  ground  state  and  hot  band  photoelectron  spectra  of  1,2,3-triazolide 
each  required  NMM  ~  109.  This  simulation  differed  from  our  previous  AFOSR  sponsored 
simulations  of  pyrrolide34  and  pyrazolide35  in  two  fundamental  aspects.  While  pyrrolide  required 
two  states  and  pyrazolide  required  three  states  in  the  simulations,  triazolide  required  four  states 
to  simulate.  In  addition  while  the  quadratic  vibronic  coupling  model  performed  admirably  for 
pyrrolide  and  pyrazolide,  the  fourth-order  quasi-diabatic  Hamiltonian,  Hd,  discussed  above,  was 
required  for  triazolide.  To  carry  out  these  calculations  three  algorithmic  extensions  were 
required,  including  the  higher  order  Hd  algorithm  described  in  section  IIB  and  the  general  open 
ended  Lanczos  algorithm  described  in  section  IIC.  We  also  introduced  into  the  Lanczos 
algorithm10  the  capability  treat  hot  band  spectra  while  reducing  the  size  of  the  vibrational  basis 
used  in  eq.  (2).  This  required  calculating  Franck-Condon  overlaps  of  the  form  (m'|n)  where 

|m')  is  hot  band  state  of  the  anion  with  one  phonon  excited,  mj=  dUj  and  n>  is  given  in  eq.  (2c). 

This  extended  our  existing  algorithm,12  which  was  restricted  to  the  vibrational  ground  state  of  the 
anion,  that  is  m  =0. 

For  the  triazolide  photoelectron  spectrum,  the  simulation  based  on  the  quartic  Hd  is  far 
superior  at  reproducing  the  main  features  of  the  experimental  spectrum,  when  compared  with 
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standard  quadratic  Hd  based  simulation  and  allowed  for  the  assignment  of  the  main  vibronic 
transitions  observed  in  the  experimental  spectrum.  Sequence  (vibrationally  excited  anion)  band 
contributions,  corresponding  to  one  quantum  in  each  of  four  low  frequency  modes  of  the 
triazolide  anion,  explain  the  observed  vibronic  transitions  occurring  to  the  red  of  the  assigned 
band  origin  of  the  neutral  species  and  some  weak  transitions  to  the  blue  of  the  origin  band.  See 
Fig.  2  where  color  codingn  is  used  to  identify  hot  bands. 

The  peak  labeled  b  occurring  at  ~3 .500-3 .525eV  (~500-621  cm"1  above  threshold  and 
~3200-3100  cm'1  below  the  minimum  energy  conical  intersection),  is  relevant  to  the  issue  of 
nonadiabatic  effects  intruding  into  what  would  be  expected  to  be  adiabatic  spectral  regions.  The 
peak  labeled  b  is  actually  composed  of  two  weak  transitions,  denoted  b'  and  b".  Such  weak 
peaks  are  frequently  assigned,  as  combination  bands  in  the  adiabatic  approximation.  However 
they  can  also  be  (nontotally  symmetric)  fundamentals  of  nonadiabatic  transitions.  In  this 
example  both  situations  are  found.  The  b'  eigenstate  arises  from  nonadiabatic  interactions  while 
the  principal  contributor  to  the  b"  level  is  an  overtone  peak.  Here  assignments  based  solely  on 
energy  can  be  problematic  owing  to  the  practical  limits  of  the  frequencies  obtained  from  ab  initio 
calculations. 

Thus,  as  in  the  case  of  pyrrolyl,  nonadiabatic  interactions  are  demonstrable  well  below 
the  energies  of  the  minimum  energy  conical  intersection.  However,  under  the  experimental 
conditions,  we  also  found  that  the  size  of  the  contributions  from  hot  bands  and  weak  adiabatic 
(allowed)  transitions  were  comparable  to  those  of  the  nonadiabatic  transitions,  complicating  the 
experimental  assignment  of  the  observed  transitions.  As  noted  above  the  existence  of 
nonadiabatic  effects  in  photoelectron  spectra  in  regions  energetically  well  below  the  minimum 
energy  conical  intersection  raised  the  question  of  the  generality  of  this  observation.  This  issue  is 
addressed  in  our  AFOSR  renewal  proposal. 

IV.  Conical  Intersections  and  Photo-induced  Decomposition  of  Energetic  Materials 
A.  Furazan  (CH)2N20 

Furazan-based,  high  nitrogen  content  energetic  materials,  especially  DAAF  (3,3 
diamino-4,4 r-azoxy furazan  =  two  amino  furazan  rings  joined  by  an  azoxy  group),  have  recently 
received  significant  attention  for  applications  such  as  explosives,  fuels,  and  propellants  due  to 
their  good  heat  resistivity,  high  heat  of  formation,  low  sensitivity,  and  good  detonation 
performance.44  The  basic  unit  of  these  materials  furazan,  pictured  above,  is  a  cyclopentadienyl 
ring  with  three  CH  groups  replaced  by  two  nitrogens  and  one  oxygen.  For  these  furazan  based 
materials  photo-induced  dissociation  of  the  electronically  excited  species  is  relevant  to  their 
denotation.  Consequently  the  photodissociation  of  DAAF  and  related  model  systems  including 
furazan  itself  have  been  studied  both  experimentally  and  computationally  by  Bernstein  and 
coworkers.44 

Our  ongoing  study  of  furazan  photodissociation45  builds  on  the  experimental  and 
computational  work  of  Bernstein  and  on  the  insights  acquired  in  our  studies  of  the  azolides. 
Bernstein's  experiments  show  that  the  initial  product  of  the  photodissociation  is  NO  and  that 
rotational  energy  in  the  NO  fragment  depends  on  the  excitation  wavelength.  Bernstein  attributed 
the  differences  in  rotational  energy  in  the  products  to  the  existence  of  distinct  pathways  through 
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different  conical  intersections.  Our  studies  of  photodissociation  in  systems  such  as  CH2COH 
+hv  ->  CH2CO+H  (Ref.  46,47)  revealed  that  g-h  plane  routing  can  also  lead  to  different  energy 
distributions  in  the  products. 

Numerical  studies  of  the  photodissociation  of  furazan  are  complicated  by  the  fact  that 
four  states  must  be  considered  in  the  entrance  channel  -  not  unlike  the  situation  in  triazolyl  - 
while  only  two,  the  states  correlating  with  NO(2II),  must  be  considered  in  the  product  channel. 
To  treat  this  properly  we  are  implementing  a  generally  weighted  state  averaged 
multiconfiguration  self  consistent  field,  GW-SA-MCSCF  procedure.48  The  need  to  consider  four 
states  in  the  reaction  channel  is  a  consequence  of  the  existence  of  two  nitrogen  atoms  in  furazan. 
This  in  turn  suggests  the  possibility  of  three  state  conical  intersections  in  the  Franck-Condon 
(FC)  region.  Indeed  in  preliminary  calculations  a  point  of  intersection  of  Si,  S2  and  S3  has  been 
located,  which  is  only  slightly  (-2000  cm'1)  higher  than  the  minimum  on  the  Si  potential  energy 
surface.  This  three  state  intersection  has  a  five  dimensional  branching  space49  which  facilitates 
motion  from  the  planar  structures  in  the  FC  region  to  nonplanar  structures. 

Using  multireference  configuration  interaction  wave  functions,  we  have  located  a  low 
energy  conical  intersection  similar  to  one  determined  by  Bernstein.  This  intersection,  Fig.  (3a), 
is  -  20,000  cm'1  below  the  energy  of  the  first  excited  state  (Si)  in  the  Franck-Condon  region  and 
lower  in  energy  than  the  azirene  (Fig.  (3c))  +  NO  asymptote  Our  ongoing  analysis  of  molecular 
motion  following  an  encounter  with  this  conical  intersection  reflects  g-h  plane  routing,  Fig.  (3a). 


g  (a)  h  (b)  (c) 

Figure  3  (a)  Two  views  of  the  furazan  low  energy  conical  intersection,  with  g  and  h  vectors  superimposed 
respectively,  (b)  an  intermediate,  and  (c)  azirene  a  decomposition  product  (CH)2N 

As  seen  from  the  g-h  plane  in  Fig.  (3a),  not  considered  in  previous  studies,  both  +g  and  ±h 
motion  produce  a  torque  on  the  NO  moiety.  We  find  motion  starting  along  -g  leads  the  local 
minimum  structure  pictured  in  Fig.  (3b).  Again  this  structure  was  not  previously  reported. 
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